Identities for density functionals linking functionals of different densities. 
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In electron density functional theory formal properties of density functionals play an 
important role in constructing and testing approximate functionals. In this paper it 
is shown that the equation 

is satisfied by a number of functionals that are of interest in density functional the- 
ory. In particular it is satisfied by E^^^ [p] , the sum of the Hartree and exchange- 
correlation energies. The parameter 7 scales the mutual electron Coulomb interaction 
energy while pn (r) and p]^_i (r) are iV-electron and (A^ — l)-electron densities de- 
termined from the same adiabatic scaled external potential of the A^-electron system 
at coupling strength 7. 
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I. INTRODUCTION 



The Kohn-Sham (KS) formulation^ of Density Functional Theory (DFT)^ is one of the 
most important tools for the calculation of electronic structure of molecules and solids. In 
all practical applications of DFT, however, approximations to the exact functionals have 
to be made^~- Exact relations for density functionals and density functional derivatives 
can play an important role in the development of accurate approximations to the exact 
functionals. A successful approach to the design of improved approximate density functionals 
is by 'constraint satisfaction'-, where the approximate functionals are required to satisfy 
properties of the exact functionals. With this in mind, the following relations for DFT 
energy functionals are derived: 

VI [pn] - VI {pI_^ 

= Jd'r{p,{r)-pl_,ir))^-J^ (2) 
[pn] - 

^/^V(p.(r)-/._,(r))^ (3) 

T2 [pn] - T2 [pI_,] 

= I dh' {p, (rO - pl^, (rO) ^igl. (4) 

r° [p^] - T° [pi_,] 

^/^V(p.(r')-p;_,(r'))|^. (5) 

ELc [pn] - El^ [pl.,] 

= I d\'{pAr')-pl.,{r'))^-§^. (6) 
J OPN irj 

Here pn and p1f_i are the ground state charge densities of an interacting and (A^ — 1) 

electron system with the same single particle multiplicative external potential fj^t ([Pa^])- 
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The potential v2y± {[pn]) is constructed to keep the charge density of the electron system 
independent of the coupling strength parameter ^Si^— that scales the electron-electron inter- 
action strength. At 7 = 1 full strength Coulomb interaction between electrons is included 
and the external potential v^y^^ {[pn]) is the external potential of the fully interacting system, 
while 7 = corresponds to the non- interacting Kohn-Sham potential. F^[p\ is the sum of 
the kinetic energy functional T"^ [p] and the mutual Coulomb interaction energy 714"^ [p] at 
interaction strength 7 and density p. [p] is the correlation part of the kinetic energy and 
^hxc {p\ i^ sum. of the Hartree and exchange-correlation energy. 

As corollaries recursion relations are satisfied. For G"' [p] = F'^[p], T^[p], V^, [p] , [p] , 
T' [p] or El^ [p] , 

M - [pl_,,] 

= d'r {pl_, (r) - p;„,„, (r)) -^^f^. M < N - 1. (7) 

Here the densities (r) are the (A^ — L)-particle densities derived from the groundstate 
wavefunctions of (A^ — L)-particle Hamiltonians with mutual Coulomb interaction 7V^e and 
external potential v]^^ {[pn]) ■ For G'^ [p] = [pn] , [pn] or El^^ [pn] , [pj] = and 
the functionals can be written as 

[pn] = J:[ dh [pl_, (r) - pl_,_, (r)) , ^^^f (8) 
L=o "Pn-l k'^) 

These identities relate functionals at different particle numbers and are likely to place 
stringent restrictions on the possible forms that potential approximate functionals can as- 
sume. 

II. PROOF OF EQUATIONS 1-3 

In the adiabatic connection approacb^^— of the constrained minimization formulation of 
density functional theory-^^ii^i^ the Hamiltonian H'^ for a system of A^ electrons is given by 

H^ = f + ^te + vl,^, [pn] ■ (9) 

Atomic units, h = e = m = 1 are used throughout. T is the kinetic energy operator, 

1 ^ 

T = ~T.'^l (10) 

i=l 
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and 'jVee is a scaled electron-electron interaction, 

N 



Kj ■> 



The the external potential 

N 

M = E^ext(M;r.), (12) 



is constructed to keep the charge density fixed at pjv (r), the ground state charge density of 
the fully interacting system (7 = 1), for all values of the coupling constant 7. The external 
potential has the Ioitdf^'^ 

vIA[Pn] ; r) = (1 - 7) Vh^{\pN];v) 

+ ^;,^(M; r) - v:{[pj,]- v)+vI,{[pn] ; r), (13) 

where vI^^^{\pn]]y) = t'ext (r) is the external potential at full coupling strength, 7 = 1, 
and 'i'g^t([PAf] '1^) is non-interacting Kohn-Sham potential. The exchange plus Hartree 
potentiaU^'^ VhxilPN]','"'), is independent of 7, while the correlation potential t'2'([pAf]; r) 
depends in the scaling parameter 7. 
The chemical potential 

= Eliv2^Ap]) - El^AvlM (14) 

depends on the asymptotic decay of the charge density^"—, and hence is independent of the 
coupling constant j^^. In Eq. iHM E'J^_-^ is the groundstate energy of the (A^ — l)-electron 
system with the same single-particle external potential f J^t {[pn] ', r) as the A^-electron sys- 
tem: 



PmI 



H^ = f + 7v;e + t)L,ext [p^ 



-^M.ext [PN\ 
M 

^M,ext [Pn] = ^^-t ( [^^] ' ^'^ (1^) 

i=l 

Note that by construction of v^^t {[Pn] ; r) , Eq- ( |T3ll . pn = Pn is independent of 7, but the 
groundstate density of the (A^ — l)-electron system derived from the same external potential 
at couphng strength 7, Pn_i, will in general be a function of 7. 
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Define the energy functional^'^^'i^ 

= 



p 



min ( 



p 



(16) 



According to tlie Levy constrained minimization formulation^'^, tlie wavefunction |^^) yields 



^ij . F'^ [p] can be decomposed as^^' 



the density p and minimizes 
The correlation energy E2 [p] is defined as^ 



1607 



(17) 



^7 
P 



- vl> 



(18) 



where |\E'°) is the Kohn-Sham independent particle groundstate wavefunction that yields 
the same density as the interacting system at coupling strength 7. Ehx [p] is the sum of the 
Hartree and exchange energy 



Ehx [p] 



and the kinetic energy functional [p] is given by. 



The full kinetic energy 



T° [p] 



T< [p] 



T 



T 



Vl/T 

P 



(19) 



(20) 



= r° [p] + t: [p] , 

with the correlation part of the kinetic energy defined as 



[P] 



p 



T 



T 



Assuming that F'^[p] is defined for non- integer electronsi^ii^i^i, at the solution point 



5pN (r) 



+ ^cxt(M;r) = p 



(21) 



(22) 



(23) 



Note that by definition of F"'[p] 

El {v2^, [pn]) = + J d'rp^ (r) t;^ {[p^] ; r) 

^^^-1 {v2.t [Pn]) = F^[pI^^] + / d\pl_^ (r) v2,t {[Pn] ; r) 



(24) 



From ([HD and 

F-'M - F^[pl^,] = f,-jd^r [pr, (r) - pl_, (r)) v^, ; r) . (25) 

Since 

I £r{p^{r)-pl_,ir)) = l, (26) 
it follows from ([23]), ([25]) and ([26D that 

F'Ipn] - F^[pl_,] = I d\ {p^ (r) - pl_, (r)) (27) 

In a recent paper, starting from the virial theorem for the interacting system^^, the author 
showed that 

Combining ([27]) and <^ yields 

[pn] - [pI.,] 
= /^^-^^U(r)-p;-Jr)). (29) 

III. PROOF OF EQUATIONS 4-6 

The mutual electron-electron repulsion energy functional can be expressed as 

VI [Pl] = E,^ [p] + ~ , (30) 

7 

which follows from Eq. ^ and the definitions ^ and fl2B^'^. Substitute Eq. ([30]) in 
Eq. ( [28]) . take the derivative with respect to 7: 

m [pn] - T2 M) - {e: [pI^,] - T2 K-i])) 

= ~- [ d^r ^^^~^ {iVhx {[pn] ; r) + {[pn] ; r)) 
7 J ^7 



7 7 97 5pAr (r) 

+ j dh' [p^ (r') - pl_, (rO) [p^] - T2 M . (31) 
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It follows from the definition of E2[pn\, Eq- (flSj) that^ 



|^i?Z[p^] = ^ {E2[pn] - . (32) 



With the aid of Eqs. (EI]), (lATOll and (El]), it can be shown that 



|dV(p.(r')-p;_,(r'))^!gl. (33) 



Integrating with respect to 7 and taking into account that T° [p] = (see Eq. (B^ ). it 
follows that 



= / d^r' [p, (rO - pl^, (rO) ^Igl. (34) 
As a consequence, from Eqs. ( |29!) and (134)) 

= /rfV(p^(r')-p;_,(rO)f^. (35) 

With 

ELc [P] = lE,. [P] + E: [p] (36) 
it follows from Eqs. ([28]), (|30]) and (|34]) that 

EL. [Pn] - El^ [pl^,] 
= J d?r' (p^ (rO - pl_, (rO) J^fLc [Pn] ■ (37) 

IV. PROOF OF RECURSION RELATIONS, EQUATION 7 

The identities in Eqs. (II]) to ([6]) all have the form 

[Pm] - [pl_,] 

= I £r' [p^ (r') - pl_, (rO) J^f' IPn] ■ (38) 



These identities were derived with 7 = 1, full interaction strength, as the starting point. The 
arguments used before are equally applicable when the Coulomb operator jVee is replaced 
by /37Ke, the coupling strength is scaled by /3 and the external potential at full coupling 
strength x = 1, is taken as v2^^{[pN]'l''^i) ■ Here /? plays the role that 7 played before. 
Recursion relations can now be derived at coupling strength 7. If the single particle external 
potential is kept fixed, i.e. for an M-electron system 

M 

v^C [Pn] = v2^t {[Pn] ; Ti) (39) 

i=l 

and pIy^ = p1.j (tij^t [Pn]) is the M-electron density constructed from a groundstate of the 
M-electron Hamiltonian H'lf^ = T + jVee + '^'^t''^ [pn] ■ Thus, from Eq. fl5S]) . but starting 
from p^_i, 

= / d\ [pl_, (r) - pl_, (r)) ^^y^} - 
Continuing this pattern leads to 

= j: ^p'^-^ - p'^---^ (^)) 1-^ M<N-i. (41) 

For G^ [pn] = VZ [pn] , [pn] or El^ [p^] , G^ [pj] = 0, hence for [p^] , [pn] and 
[pn] 

L=0 ^Pn-L l^J 

V. DISCUSSION AND SUMMARY 

The identities reported here were derived by reference to the eigenf unctions of many 
particle Hamiltonians, though indirectly, via Eq. f l23p . They therefore are valid for w- 
representable densities^*^, in other words for densities that can be determined from ground- 
state wavefunctions of a many-particle Hamiltonian. The assumption was made that all 
functional derivatives are well behaved and this implies that the functionals are defined for 
non-integer particle numbers^l. In the derivations no explicitly reference was made to the 
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structure of the wavefunctions or density, the equations are therefore vahd for all groundstate 
densities. Specifically they are valid for all degenerate densities. 

Equation may be the most useful of the set presented here for testing approximate 
exchange correlation functionals. The obvious test can be done at full coupling strength, 
where for a real system the external potential is the Coulomb interaction with the nuclei, 
and requires two separate self-consistent Kohn-Sham calculations for the and (A^ — 1) 
electron systems to determine pj^ (r) and p]v-i (^) • -^^^ [p\ ^ single Kohn-Sham calculation 
is sufficient. The other identities are more difficult to use for testing purposes since they are 
satisfied by the functionals of the exact densities which will require independent calculations 
to determine accurate densities. 

If the functionals can be expanded as a series in the usual form^^ 

[p] 

= G^[po] + f^^ [ rfVi...dV„(p(rO-po(ri))...(p(r„)-po(r„.)) ^ ("T [^'] ^ ^^^^ 
^n!J dp{ri) ...dp{rn) 

then ( l38ll implies that 

»=f:^J ^n.^-'^r-' (^.) - (^0) ..^ (P. M - (r„)) .JIZ^Iky 

(44) 

for F'^'Ip], T'^Ip], V^, [p] , T^' [p] , T° [p] and El^^ [p] . Equation ( l44l) in itself places a constraint 
on the structure of potential approximate density functionals. 
Note that for 7 = 0, 

P%-L{r)-p%-,-L{r) = \<PN-L{r)\' (45) 

where (pi (r) is an eigenfunction of the Kohn-Sham Hamiltonian with corresponding eigen- 
value Ei ordered so that < ej+i. Equation ( H2l) . for 7 = 0, becomes 

[pn] = e y \^^-^ — S^^Tw — ' ^^^^ 

and p%_L (r) = ^i^"' 10. (r)l' • Eq. m is valid for G^' [p^v] = ^Lc [Pn] and for G^ [p^] = 
^^hxc [Pn] and since lim^^o^-E'2 [pn] = 0^, it follows from ( l36l) that 

(47) 
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N-2 

Ehx [pn] = 5Z / 

L=0 ^ 



d^r\<j)N^LM\ Vm 



N-L 

E 

i=l 



From the definition of tlie Hartree potential, 




/ 1 ' 
r — r 



it follows form Eq. (|6]) and fl36|) that 



Ex [pn] — Ex [p%-i\ 

j dhvx{\pN]:r)\(l)N (r)l' 

+1 / dhdh'\<t>N{r)?-^MN{rf. 
2 J |r — r I 



(48) 



Eqs. fj47j) and f l48|) were discovered by Levy and Gorling some time ago.— 

The only other identity that survives at 7 = is the independent particle Kohn-Sham 
kinetic energy: 



This expression follows trivially form the Kohn-Sham expression for the independent particle 
kinetic energy. All identities are valid at 7 7^ 0. 

In summary, equations that a set of exact density functionals satisfy for functionals of 
different densities, where the densities are derived from the same extern potential, were 
derived. As a corollary, it was shown that the functionals can be expressed as a sum over 
integrals of functional derivatives, where the sum runs over all if-representable densities for 
the same external potential that integrate to an integer less than the particle number. These 
relations place stringent constraints on functionals that appear in density functional theory 
and it will be difficult to satisfy by approximate functionals. 

Appendix A 



T° [pn] - T' 





(49) 



FromEq. ([HD 



E2 K-i] - 



K-i] 



7 




(Al) 
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The derivative of E] [plr_i\ with respect to 7, from definition ([T8|) . can therefore be expressed 



as 



^^^^ [^^-i] 
7 



7 



(^7 /'iv-i 



T + 7V; 

T + 7V; 



7 

Pn-1 



9 T7 

— 

Pn-1 



Pn-1 



d_ 

' Pn- 



.^0 



(A2) 



Upon adding and subtracting (c.c. stands for the complex conjugate of the previous term) 

d 



+ 



(^7 Pn-1 
(^7 Pn-1 



Pn- 



+ C.C 



7=0 



K-l,o.t[pN]\^\_^ )+C.C 



and utihzing the normahzation of the wavefunctions which imphes that 



d 



(^7 \ ''JV-I ''JV 



(A3) 



(A4) 



Eq. fl22l) becomes 
d_ 



El K-i] 



E2 K-i] - T2 K-i] 



7 



^7 



(^7 P]v- 
d 



' Pn-1 



Pn-1 



57 \ ^ 

This equation can be simphfied since 
d 



Pn-1 



(A5) 



Pn-1 



Pn-1 



\v2^t [P- 



•N 



9 T7 

^7 ''at-i 



^7 



^^ext (M ;r) 



(A6) 
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^0 



Ehx [pn-i] , 



(A7) 



where E^x [p]^_i\ the sum of the exchange [Pat-i] and Hartree interaction energy 
U [p'li_i\ of the (A^ — l)-electron system and the charge density pjj^i is a function of 



7^, 



where 



"±'7 

Pn-1 



(Tr 



dpl-i (r) 
d'y 



Vhx{[pN.i\ ;r) , 



•^P]^-! (r) 



(A8) 



(A9) 



) and 



y yield the same density p^_i 



is the sum of the exchange and Hartree potentials for the (A^ — l)-electron system. Using 
Eqs. (jMl), (M, (USD and the fact that 
by construction, Eq. flASP can be expressed as 

^^^^ K-J 

_ E2 K-i] - K-i] 



7 



(97 



(^J ([Piv] ; r) + -fVhx {[Pn] ; r) - -fVhx {[Pn-i] ; r)) (^^0) 
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